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We consider Gaussian signals, i.e. random functions u(t) (t/L £ [0, 1]) with independent Gaussian 
Fourier modes of variance ~ l/g a , and compute their statistical properties in small windows [x, x+5]. 
We determine moments of the probability distribution of the mean square width of u(t) in powers of 
the window size 5. We show that the moments, in the small-window limit S <C 1, become universal, 
whereas they strongly depend on the boundary conditions of u(t) for larger S. For a > 3, the 
probability distribution is computed in the small-window limit and shown to be independent of a. 



Gaussian signals — random functions with independent 
Gaussian Fourier components of variance ~ l/g" — have 
been used to describe physical situations ranging from 
1/f noise in electric circuits Q to intermittency in tur- 
bulent flows , and to interfaces in random media 0, • 

The exponent a > 1 of the power spectrum fixes the 
signal's mean square width W2(L) on a scale L, as char- 
acterized by the roughness exponent £ 
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£ 2C , (1) 



which, in one dimension, is given by C = Ua - 1). The 
notorious random walk, with roughness £ = i, corre- 
sponds to a power spectrum with a — 2. The well-studied 
curvature- driven model Q also belongs to this class of 
systems and corresponds to the case a = 4. 

For elastic interfaces in disordered media and other 
systems with non-integer cv, this exponent of the power 
spectrum and equivalently the roughness exponent are 
extremely difficult to calculate 0, U@, EJ ■ Exact results 
are rare E3- I n contrast to the random walk and the 
curvature-driven model, these systems are not exactly 
Gaussian. This means that the Fourier modes, even at 
small q, are correlated. More intricate statistical prop- 
erties than the scaling of u>2 are able to expose these 
correlations. A common approach consists in studying 
the probability distribution of the mean squared width 
W2 , which fluctuates from sample to sample 0, flil ITsT ] . 
This distribution has been used to characterize the ge- 
ometric properties of numerical and experimental data 
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Non-Gaussian corrections for the probability distribu- 
tion of the mean-square width were explicitly determined 
in a non-trivial depinning problem and found to be on the 
0.1% level 13]. It was possible to understand this small 
effect because non-Gaussian corrections appear only in 
high orders of perturbation theory 

[HEil- The excel- 
lent agreement between complicated physical models on 
the one side and their effective Gaussian description on 
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FIG. 1: Periodic (left) and free (right) Gaussian signals 
u(t-L). This paper studies the statistics of signals inside 
a small window (middle). 



the other motivates a finer analysis of the universal sta- 
tistical properties of Gaussian signals, which is the object 
of this letter. 

Because of their definition in Fourier space, it is most 
convenient to study periodic signals with u(t) = u(t + L). 
However, experimental systems are usually non-periodic. 
Free boundary conditions arc commonly modeled by 
Gaussian signals u(t) with zero mean and vanishing 
derivatives at the end-points (see Fig. The proba- 
bility distribution of u>2 depends on the boundary con- 
ditions: The different distributions of free and periodic 
signals were computed analytically 0, 0] . 

Several authors have studied the signal u(t) inside a 
small window (see Fig.^), i.e. the piece with t/L g [x, x+ 
5]. Antal et al. |l9l| performed numerical simulation for 
different values of a. It was empirically found that the 
probability distribution of the mean squared width inside 
a window agrees well with free boundary conditions for 
a ~ 2 15]. The two distribution differ markedly for a 
outside this range, as was clearly shown by de Queiroz 

In this paper, we calculate analytically the statisti- 
cal properties of Gaussian signals, with periodic or free 
boundary conditions, in small windows: We determine 
the expansion in powers of the window size S of the lower 
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moments of the probability distribution of W2 and find 
them to be identical for the periodic and for the free sig- 
nal far from the boundary. This is important for the 
analysis of experimental data, for which the boundary 
conditions can usually not be controlled: the statistics 
in a small window is independent of them, i.e. contains 
only universal information. We complement this deter- 
mination of moments by numerical calculations of the 
probability distribution itself. 

The Gaussian signal is defined by the action 



b ~ 2 



dt 
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where the derivative for non-integer a is understood in 
Fourier space [l9ll2l|. For a periodic signal of zero mean, 
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this means that the a n and b n are Gaussian random num- 
bers of variance c 2 = L 2< "2 _2 '» • (7m)~ Q . Free Gaussian 
signals are commonly modeled by cosines with period 2L, 



U(t) = ^ °n COS {^-t 

n=l 



(4) 



where the c„ are of variance 2L 2 ^ • (irn)~ a . 

The mean squared width of a signal u(t-L) in a window 
[x, x + S] is 



w 2 (x,S) = - 
o 



-S / 1 f x+8 \ 

u 2 (t-L)dt- I- J u(t-L)dt \ 
For the free boundary conditions, we have 

Wzix, d) = ^2 c n c m D nm (x,S) (5) 

with 
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D nm (x,S) = — / cos(7rnt) cos(7rmt)drj 
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The periodic signal's mean squared width, Wj™ > has an 
analogous expression in terms of cos-cos, cos-sin and sin- 
sin integrals. Clearly, w^ 1 does not depend on the initial 
point x, but only on the window size 6. 

The above equations allow to compute W2 for one given 
sample (choice of {a ni b n } or {c„}) and its probability 
distribution, which is characterized by the ensemble av- 
erage (u>2) and by the rescaled distribution <j>(z) with 



z = W2/ (^2). For free boundary conditions, Eq. (JSj im- 
plies 



2£ 2C y> Dnn(x,6) 



This gives for the variance of the rescaled probability 
distribution (f>(z) 
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which is independent of L. All the other cumulants of 
4>(z) are scale-free, and are defined analogously through 
multiple sums. We note that for 5 = 1, and with the 
normalization l|2|l. we have 

= i)> = ™C(oO 
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where ((a) is the Riemann zeta function. All the (scale- 
free) cumulants are known for 6=1 0, 0] 
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Sums as in Equations © and J7J) may be evaluated 
with a powerful formula [2l| 
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where [a\ is the integer part of a. The Eq. Qlljp. which 
is in the spirit of the Euler-Maclaurin formula, is valid 
only for analytic functions / and non- integer a: the first 
term on the right can be interpreted as the naive limit of 
the sum as 5 — > 0, with t = nS. The second term on the 
right contains the Taylor expansion of f(nS) around zero. 
For integer a, the singularity of £(1) generates additional 
logarithms, which are explicitly known |2*lj . 

We use Eq. 1)10(1 . and its generalization for integer a, 
to compute moments of the probability distribution of 
w 2 in a window [x, x + S] (for a periodic signal, the result 
is independent of x, for a free one, this is true only for 
a = 2). For periodic boundary conditions, we find 
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This formula gets modified by logarithms for a = 3, 
a = 5, etc |2l|. Interestingly, an analogous expansion 
appears in the correlation function governing the density 
of zero-crossings of a Gaussian signal j2^. We note that, 
for even a, the series (fTTf) stops, because the Riemann 
zeta function vanishes for even negative integers. For ex- 
ample, the periodic random walk (a = 2) and the driven 
curvature model (a = 4) yield 
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From the series (|11[) . the dominant term of (w^ 
scales for small windows as (L(5) a_1 = (L<5) 2< >, in agree- 
ment with the self-affinity relation (JTJ. However, this 
natural scaling breaks down for a > 3, i.e. for rough- 
ness exponents £ > 1. There, the small window scaling, 
from Eq. ljTT|) . is as (LS) 2 ■ L 2 ^ 2 . In addition to an 
a- independent scaling at small distances, a scale factor, 
depending on a and on the system size appears explicitly. 
This was pointed out by Tang and Leschhorn |23j in the 
context of depinning, where models with £ > 1 (a > 3) 
appear naturally [24 l2fl|. 

For a < 3, and free boundary conditions with < x < 
1, the dominant integral in Eq. I|1U|) involves an intri- 
cate double limit, where (5 — > and x/S — ► oo. Rapidly 
oscillating terms with vanishing contributions, generated 
from the second limit, need to be discarded. We find [2l| 
that the dominant term of (w^) proportional to S a ~ 1 is 
identical to the dominant term of (w^ 1 ), from Eq. (jllfl . 
The mean squared width (1U2) is thus insensitive to the 
boundary conditions. We note that the expansion of (w^) 
obtained from Eq. (jlOJl provides non-intuitive explicit 
prescriptions for extracting the roughness exponent from 
experimental or numerical data in powers of 5, which dif- 
fer from the standard ansatz [2^. 

For a > 3, the dominant term for small 5 is 
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The vanishing derivatives at the end-points force (iz4 r ) 
to be smaller than (w^ 01 ) for all < x < 1. At the 
end-points x — and x = 1, the <5 2 term vanishes [211 . 

We now turn to the higher moments of the distribution 
4>(z), (cf. Eq. J7J) for free boundary conditions), which 
can be evaluated with straightforward generalizations of 
Eq. (1KH to multiple sums. The leading term of K n , for 
a < 3, is given by an n-dimensional, well-behaved, in- 
tegral, which can be computed numerically. The calcu- 
lation of the variance K2 is along the lines of the above 
determination of (u^)- Discarding again rapidly oscillat- 
ing terms in the limit S — > and x/8 — * 00, we can an- 
alytically prove that k^ c1 (8) coincides with k,2(x,8) 
for 1 < a < 3. Furthermore, in the special case a = 2, 
the integral explicitly agrees with the variance K 2 r (8 = 1) 
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FIG. 2: Probability distribution 4>(z) for free and periodic 
signals at a = 2.5. Direct simulations for 8 = 0.01 (with 
x = 0.495 for the free case) are compared to analytic solutions 
for 8= 1 EH]]!. 



from Eq. which for the random walk is independent 
of 8. This is due to the Markov-chain property of the free 
random walk. The density matrix generating the path is 
infinitely divisible, and essentially reproduces itself for all 
values of 8. In contrast, for all a ^ 2, we have 



0) = 0) ^ = 1). 



(14) 



'Windows' of periodic or free signal are thus different 
from the free signal itself. To stress this crucial result, 
we show in Fig. [21 the numerically computed distribution 
<p(z) for 8 = 0.01 for both free and periodic boundary 
conditions, and compare them to the distributions for 
8 = 1. The two distribution at 8 — 0.01 agree very well, 
but, for this value of ( close to |, they differ somewhat 
from the free 8=1 distribution. This point was missed 
in previous work [l^, because the exponent considered 
there was very close to i. Notice that the distribution 
for the periodic 8 = 1 signal is very different from the 
other distributions. 

The calculation, for a < 3, of cumulants n n with 
n > 2 presents no conceptual difficulties, and we conjec- 
ture that all moments, and thus the distribution itself, 
in the small- window limit, are independent of boundary 
conditions. This is also supported by the numerical re- 
sults of Fig. □ 

In the case a > 3, the integral term in Eq. (|10|> is sub- 
dominant and the calculation of all cumulants becomes 
feasible. We find for the moments 
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This implies that the distribution 
window limit is explicitly given by 
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FIG. 3: Probability distribution <f>(z) 
3.5. A direct simulation of Eq. (jj) (i 



for free signals at a = 
= 0.45, 5 = 0.1) with 



10 6 samples is compared to the function I16H in the small- 
window limit and to the analytic solution for 5 = 1 [Iflj . 



The exact small- window distribution of the mean-square 
distribution, valid for all a > 3, was already known to 
be exact in the limit a — > oo . In Fig. El the small- 
window distribution is compared to the numerically ob- 
tained histogram for a — 3.5. Notice that, for a > 3, a 
small window of the free signal has completely different 
statistics from the free signal itself. 

In conclusion, we considered in this paper statistical 
properties of Gaussian signals. We studied the influence 
of boundary conditions on the signal in a small window. 
An exact sum formula, non-trivial generalization of the 
Euler-Maclaurin equation, allowed to reach a systematic 
procedure for computing moments of the mean square 
distribution function which, at small S, was found to be 
independent of the boundary condition. 

The expansion in powers of the window size 6 changes 
at a — 3 (corresponding to a roughness of £ = 1). Above 
this value, the calculation of all moments of the distribu- 
tion function (in the 8 — > limit) becomes particularly 
simple, and the whole probability distribution was com- 
puted. Clearly, the small-window limit studied in this 
paper plays an important role: It is independent of the 
boundary conditions and contains the true universal in- 
formation a Gaussian signal. 
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